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SLENDER BODY TREATMENT OF SOME SPECIALIZED PROBLEMS ASSOCIATED 


WITH ELLIPTIC-CROSS-SECTION MISSILE CONFIGURATIONS 

Raymond L. Barger 
Langley Research Center 


SUMMARY 

Slender body methods have been applied to some specialized problems 
associated with missile configurations with elliptic cross sections. Expres- 
sions are derived for computing the velocity distribution on the nose section 
when the ellipse eccentricity is varying longitudinally on the missile. The 
cross-flow velocity on a triform fin section is also studied. 


INTRODUCTION 

In recent years, interest has grown in missiles with bodies having ellip- 
tic cross sections. One reason for this increased interest is the possibility 
of an aerodynamically efficient mounting of such a low-profile missile on 
the parent airplane. Another advantage of such missiles over conventional 
circular-cross-section configurations is the increase in directional stability 
that occurs with increasing ellipse eccentricity (ref. 1). 

A primary tool for missile analysis is slender body theory (ref. 2). One 
of the advantages of this theory is that it permits the nose and finned sec- 
tions of the configuration to be treated separately. Thus a missile with hori- 
zontal fins, constant-eccentricity elliptic cross sections in the nose section, 
and an elliptic cylinder body in the finned section could be analyzed by com- 
bining the results of reference 2 (sec. 4.3), 3, or 4 for the nose section with 
those of reference 5 f or the finned section. 

This report presents some further results for missiles with elliptic- 
cross-section bodies. First, velocity expressions are obtained for the nose 
section with ellipse eccentricity varying longitudinally along the missile 
(fig. 1). This result is important because the shape of the equipment that 
must be carried in the nose section closely controls the cross-section shape 
and, consequently, the permissible eccentricity. The analysis is based on a 
result derived in reference 4. 

Second, expressions are obtained for the cross-flow velocities in the 
finned section for a triform configuration. This analysis utilizes, first, 
a conformal transformation of the ellipse to a circle, and then a series of 
transf ormations that has been previously used by Bryson (ref. 6). 

The present work represents a supplement to the already substantial body 
of slender body results applicable to missile conf igurations . The use of the 
expressions obtained herein, together with those already available, permits the 


analysis of certain configurations not treated in previous work. A separate 
computer program has not been written for the calculation of these expressions. 
Rather, they are to be incorporated into already existing slender body programs 
to permit application of those programs to a larger class of missiles. 
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SYMBOLS 

ellipse semimajor axis 
ellipse semiminor axis 
defined by equation (9) 

= s/a2 - B2 

see equation (35) 

see equation (39) 

y-coordinate of horizontal fin tip 
see equation (3*0 
complex force 
lift force 

free-stream Mach number 

free-stream dynamic pressure 

cross-sectional area of missile body 

coordinate directed along missile axis 

coordinate in direction of undisturbed free stream 

see equation (37) 

z-coordinate of vertical fin tip 

see equation (36) 

perturbation velocity components in s-, y-, and z-direction, 
respectively 

complex velocity potential 

complex variable in plane perpendicular to missile axis, y + iz 
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x c complex coordinate of cross-section center of area for force 

calculation, -ias 

X'j,X 2 ,X 3 ,X 4 complex variables in successive transform planes 

y,z real variables in missile cross plane 

a angle of attack 

3 = si m2 - 1 

3 S sideslip angle 

5 see equation (33) 

r|,£ elliptic coordinates (eq. (la)) 

Hq value of n representing missile section contour 

£> = n + i? 

K c defined by equation (27) 

p defined by equation (2) 

Po = P(n 0 ) 

$ velocity potential, real part of W 

Subscripts: 

h associated with side flow 

t associated with thickness 

a associated with angle of attack (vertical component of free stream 

in cross plane) 


ANALYSIS OF FLOW ON MISSILE NOSE SECTION WITH 
VARIABLE-ECCENTRICITY ELLIPTIC CROSS SECTIONS 

In this section, expressions are derived for computing the velocity com- 
ponents due both to thickness and to angle of attack on a nose section having 
elliptic cross sections that vary in eccentricity longitudinally along the 
missile axis. It is also shown that, within the slender body approximation, 
the total force on the nose section depends only on the cross-section shape at 
the base of the nose section. 
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Basic Relations 


The coordinate system is taken along the missile axis (s-axis) and perpen- 
dicular to it (x- and y-axes) (fig. 2). The velocity distribution on the mis- 
sile is more conveniently obtained with the use of elliptic coordinates n,C 
defined by the transformation (see ref. 7, pp . 157-159), 


y + iz = x = c cosh 


c cosh (n + i£) 


or 


y = c cosh r) sin 


z = c sinh 


n sin T 


(la) 


(1b) 


From equation (la), one obtains the relation, 



c sinh £ 


Another useful variable is defined by 


P 


ce 1 "! 

~2 


On the surface o = Bo , 

A = c cosh rig 
B = c sinh no 

and p has the value 


ce n ° A + B 


PO = 


( 2 ) 


(3) 


( 4 ) 


Thus the ellipse n = Ho maps into the circle p = pq. Then the transforma- 
tion (eqs. (1b)) can be written in terms of p and Z, as follows: 
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1 / \ 
y = — ce^ + cos £ 

2 \ ceV 


1 / \ 

z = — Ice 1 ^ - sin XL, 

2 \ ce^/ 


f ^ r 

p + — cos C 

^ W 


p sin C 

V 4p/ ^ 


( 5 ) 


The surface values of dp/ds and d£/ds will be needed in the analysis. To 
obtain them, first differentiate equations (5) with respect to s: 


1 - 


1 + 


a 

a 


dp 

( c 2 \ 

d? 

c 

dc 

COS XL — - 

P + — 

sin £ — = 

- — cos 

? — 

ds 

\ 4p/ 

ds 

2p 

ds 


dp 

1 c2 \ 

d£ 

C 

dc 

— + 


COS XL — = 

= — sin 

e — 

ds 

\ W 

ds 

2p 

ds 


( 6 ) 




The coefficients simplify considerably on the missile surface because of the 
relations 


1 


c2 

^Po 2 


2B 

A + B 


1 



2A 

A + B 


PO 



A 


PO 


C 2 

4 P0 


= B 
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Making these substitutions in equations (6) and solving for the derivatives 
yields, on the missile surface, 


dC 


sin £ cos £ 


d( A 2 - 

B 2 ) 

ds 

2( A 2 

sin 2 C + B 2 cos 2 

o 

ds 


dp 

A ; 

sin 2 5 - B cos 2 £ 


d( A 2 - 

B 2 ) 

ds 

4 (A 2 

sin 2 5 + B 2 cos 2 

o 

ds 



( 7 ) 


Velocity Components on Nose Section Due to Thickness 

According to slender body theory, the velocity components due to thickness 
can be determined separately from those due to angle of attack. The scalar 
potential function representing the effect of thickness is derived in refer- 
ence 4 (eq. (23)) to be (in the notation of this report) 


<Mn,c) 


1 dS(s) / c\ 

(ri 4- In - 

2 IT ds \ 2/ 


b 0 


A 2 d /B\ 
4 ds \A/ 


e -2(n-n 0 ) 


cos 2c; 


( 8 ) 


where 


1 dS(s) 8 f S 

bn = — In - - \ In (s 

2tt ds 2 J n 


d 2 S(a) 

a) da 

da 2 


(9) 


and unit undisturbed free-stream velocity is assumed. From equation (2), 


c 

n + In — = In p 
2 


and on the surface 


^ e 2n 0 

4 ds\A/ 


c 2 cosh 2 no 2 rio 
4 6 


— (tanh no) 
ds 


c 2 

r 


e 2 no 


dno 

ds 


PO 2 


d no 

ds 
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Thus when the eccentricity is constant, this quantity is zero. With these 
substitutions, equation (8) becomes 


i ds(s) , d no _ , % 

<t>f(p,5) = — In p + b 0 + p 0 2 e -2r l cos 2 % (10) 

2tt ds ds 


Then the axial perturbation velocity is determined on the surface by dif- 

ferentiating equation (10): 


• d<(> 
u t = — 
ds 


1 d 2 S(s) 1 dS(s) 1 dp 

In Pg + — 


2lT ds 2 


2tt ds Pq ds 


db 0 dp 0 


-2q 0 


dno 


n=no 


+ + 2pg e u cos 2£ 

ds ds ds 


? d2n ° -2n n P -2n n dn ° d H 

+ Pq 2 e 0 cos 2£ - 2pg 2 e 0 cos 2£ — 

d s 2 ds ds 


n=no 


o -2nn dX] ° dC 

2pQ 2 e ^ sin 2 £ — 

ds ds 


n=no 


where from equation (2), 


dn 

ds 


1 dp 

n=no ds 


n=no 


1 dc 
c ds 


( 1 1 : 


and from equation (9), 


dbg 

1 

d 2 S(s) 3 

1 ^ 

d 2 S(0) 

ds 

2tt 

ds 2 2 

ds 2 



(s 


a) 


d3s(a) 

da 

da3 
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Combining terms and substituting 


e -2n 0 


A - B 

yields 

A + B 


1 d 2 S(s) $P0 

ut = — In 

2tt ds 2 2 


1 dS(s) 1 dp 1 dS 2 (0) 

2v ds Pq ds n=no ^ ds 2 



(s 


O) 


d3s(a) 

da 3 


da + (A - B) 


COS 


dpo dr io c 2 d2 no 

2C + — cos 2r 

ds ds 4 ds 2 


c 2 

- — cos 


2 


2C 


dOo dr) 
ds ds 


n=no 


C 2 dr 10 

- — sin 2£ 

2 ds 


d? 


ds 


n=no 


( 12 ) 


where dn/ds is given by equation (11) with dp/ds and d£/ds computed from 
equations (7). 


For the special case of constant ellipse eccentricity, the last term in 
equation (8) is zero. Consequently, only- the first four terms of equation (12) 
are retained for this case which is treated, for example, in reference 3* 


The calculation of the cross-flow velocities resulting from the thickness 
effect is simplified somewhat by the use of the complex potential W^. A for- 
mula for can be obtained from equation (10) for the scalar potential by 

the method of reference 7 (p. 126). The result is 


1 dS(s) 

w t (5) = - - — 

2TT ds 


c\ 9 dr)0 or 

E, + In — + bo + Pep e - ^ 

2/ ds 


Then v^ and w^ are determined by 


Here 


v t - iw t 


dW t dC 
d£ dx 


d£ 1 1c sinh E, 

dx dx/dC c sinh t c s j_ n h \ 


c sinh n cos X, - ic cosh n sin €, 
c 2 (sinh 2 r) + sin 2 £) 


(13) 


(14) 


(15) 
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where the bar indicates complex conjugate. On the surface, equations (3) 
apply. By making these substitutions and replacing c 2 by A 2 - B 2 , one 
obtains 




B cos C, - iA sin C, 


dx 


n=no 


A 2 sin 2 C, + B 2 cos 2 £ 


( 16 ) 


The velocity in the £ -plane is 


dW t 1 dS(s) 0 dr 10 

r 2p 0 2 e" 2 ? 

d£ 2 tt ds ds 


On the surface 


dW t 

dT 


n=no 


1 dS(s) . _ 2 n n dr io 

— - 2po 2 e 0 (cos 25 - i sin 2C , ) 

2it ds ds 


(17) 


Substituting equations (16) and (17) in equation (14) yields 


dW t 

dx 


n=no 


1 dS(s) 


drio 


2po^e” r|( - ) (cos 2c, - i sin 2C,) 

2 it ds ds 


B cos C, iA sin C, 

A 2 sin 2 5 + B 2 cos 2 C, 


Thus, 


B cos C, dS(s) 
2tt ds 


c 2 dl "l0 

— (B cos 2c, cos C, - P sin 2C, sin £) 

2 ds 


v t 


n=no 


A 2 sin 2 c, + B 2 cos 2 C, 


(18) 


Wt n=no 


A sin C, dS(s) c 2 dr l0 

- — (B sin 2c, cos C, + A cos 2c, sin C.) 

2tt ds 2 ds 


A 2 sin 2 C, + B 2 cos 2 c, 


(19) 


For the special case of constant ellipse eccentricity, only the first term 
in the numerator remains for both components . 
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Velocity Components on Nose Section Due to Angle of Attack 

The perturbation velocity potential for the cross flow due to angle of 
attack is (see ref. 7, p. 159), in the notation of this report, 


W a = -ia(A + B) sinh (£ - n 0 ) 


(20) 


The scalar potential, which is the real part of W, is 




= 0t sin £ 



and its derivative with respective to s is 

3<J> d P0 3(J) d£ 3cj) dp 
ds 9 pq ds. 3? ds 9p ds 

Thus , 


2Po 

d P0 

( PO 2 ' 

\ d ? 

( PO 2 ^ 

(dp 

sin £ 

+ cos £ 

p + — 

— + sin £ 

1 1 ~ 


P 

ds 

\ p / 

'ds 

\ P 2 / 

1 ds 


On the surface, 


P = P0 

d Po 1 /dA dB\ 

ds 2 (ds ds/ 

and d£/ds and dp/ds are obtained from equations (7). Substituting these 
quantities into equation (21) and simplifying yields 


% = 


a sin L, 



dA 

— + B(A + B) 
ds 


cos^ 


dA _ 

£ — + A(A sin^ 
ds 


C 


a 2 sin^ 5 + b2 cos2 £ 


- B cos2 



( 22 ) 
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For a body whose cross sections have the same eccentricity, 


dB B dA 

ds A ds 

and equation (22) reduces to 


aA(A + B) sin £ — 
ds 

a A 2 sin 2 £ + B 2 cos 2 £ 

which is the known result for this case. (See ref. 2, p. 78, eq. (4-43).) 

The cross-flow velocity components due to angle of attack are obtained 
from the relation 


v a - iw a 


dx 


d £ 

d£ dx 


d£ 

-ia(A + B) cosh (£ - pg) — 

dx 


Evaluating this expression on the surface with the use of equation (15) yields 


v a - lw a = 


-ia(A + B) jjiosh (p - pg) cos £ - i sinh (p - pg) sin £j (B cos £ - iA sin £) 


A 2 sin 2 £ + B 2 cos 2 £ 


v a - lw a| n - 


n=no 


-ia(A + B) cos £ (B cos £ - iA sin £) 
A 2 sin 2 £ + B 2 cos 2 £ 


Thus, 


-aA(A + B) sin £ cos £ 
A 2 sin 2 £ + B 2 cos 2 £ 


( 23 ) 
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w a 


aB(A + B) cos 2 £ 

A 2 sin 2 £ + B 2 cos 2 £ 


(2H) 


These, of course, are perturbation velocity components. The total vertical 
velocity due to angle of attack is 


a 



B(A + B) cos 2 £ 
sin 2 £ + B 2 cos 2 


C 


1 


Total Force on Nose Section 

For a missile at angle of attack, the complex force from 0 to s is given 
in reference 2 (eq. (3-62), p. 50) and, in the notation of the present report, 
is 

- = 4-rrai + 2 ^=[s(S) x c (i)l (25) 

q ds u u 


where a-) is the coefficient of x -1 in the Laurent expansion of W(x) and 
where the coordinates are referred to the s-axis parallel to the free-stream 
direction. Thus, in equations (13) and (20) for the complex thickness and lift 
(angle-of-attack) potentials, x is replaced by x - x c , where 


x c = -ias 


Then in the equation for the total potential, 


W = W a + w fc 


the thickness component is written 


W t = bg(s) + — 


1 dS(s) x - x c + J(x - x c ) 2 - c 2 dn 0 ?r 

In — — + po 2 e 

2tt ds 2 ds 


( 26 ) 


where 


e~^° a — jx - x c - \j(x - x c ) 2 - c 2 
c 


( 27 ) 
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Similarly, with the appropriate substitutions,- equation (20) for the com- 
plex potential due to lift can be written 

W a = - j-jx - x c + \J(x - x c ) 2 - c 2 - [* - x c ~ \J( X ~ x c) 2 “ c 2 jj (28) 

The contribution of the first term of (eq. (26)) to a-) is zero. 

x c dS(s) 

The contribution of the second term is known to be - — . To find the 

2tt ds 

contribution of the third term, which accounts for the variable eccentricity 

_ 2 £ 

of the cross section, one must expand e 0 as follows: 


5 -2 ^ c = — (x - x c ) 2 - 2(x - x c ) \J (x - x c ) 2 - c 2 + (x - x c ) 2 - c 2 ] 


2(x-x 0 ) 2 2(x - x G ) 


(x - x 0 ) - 


2 ( x - X C ) 8(x - x 0 )3 


- 1 


0 + 


4(x - x c ) 2 


_ o r 

Thus there can be no x~1 term in the Laurent expansion of e c . Conse- 
quently the contribution of the last term in equation (26) to is zero, and 

the effect of the variable eccentricity on the total complex force is zero. 

By a similar expansion of the right side of equation (28), the contribu- 
tion of W a to a-; is found to be 


a 1 ,a 


iaA(A + B) 
2 


Summing all the contributions to a^ and substituting in equation (25) 
yields the force, which in this case is lift, 
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- = 4TT(ai t + a! «) + 2 it 

q 


d(AB) 

ds 


dx r 


+ AB 


ds 


= -2irx r 


d( AB) 
ds 


+ 2‘/TiaA(A + B) + 2'irx r 


d( AB) 
ds 


dx r 


+ 2itAB 


ds 


Substituting dx G /ds = -ia into this expression yields for the lift 
— = 2ttcxA 2 

q 

CROSS-FLOW VELOCITY IN FINNED ’ SECTION FOR TRIFORM CONFIGURATION 

If the missile has only a horizontal set of fins, the velocity distribu- 
tion in the region of the fins can be computed by the method of reference 5. 

If the missile also possesses a vertical set of fins of equal length, the 
vertical onset cross flow is not affected by the vertical fins because of the 
symmetry of the configuration. Similarly, the streamlines in the presence of 
a horizontal onset flow are not affected by the presence of the horizontal 
fins. Consequently , the effects of this horizontal component of the onset flow 
can also be computed by the method of reference 5 by switching the appropriate 
variables in that solution, in a manner somewhat similar to that of reference 8 
or reference 2 (p. 123). 

A third case of considerable interest is a triform configuration with a 
single vertical fin. For vertical cross flow (fig. 3), the potential is the 
same as if the configuration had only horizontal fins because of the symmetry 
of the flow. For horizontal cross flow, on the other hand, the flow is not 
symmetric. For this case the potential function can be found by first mapping 
the elliptical cross section into a circle and subsequently mapping the result- 
ing circle trifin configuration into a flat plate. The latter mapping can be 
accomplished by a series of transformations used by Bryson (ref. 6). 

The transformations and the corresponding contours are outlined in fig- 
ure 4. The mapping 


x 2 = -(x-, ± 7 X 1 2 - c 2 ) 


(29) 


maps the ellipse into a circle. In terms of elliptic coordinates E, □ q + i£, 


x 2 


= le? 


(30) 
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The circle is flattened by the transformation, 


x 3 


= x 2 


+ 


Po £ 


x 2 


(A + B) cosh (K - Hq) 


(3D 


where Pq is defined by equation (4). Finally the resulting perpendicular 
segment is mapped into a single vertical segment by the transformation, 


where 


x 4 


= ± ^xg 2 - d 2 - i6 


6 = — \\[t 2 + d 2 - d 
2 x 


PO^ 

d = D + 

D 


= 1 (d + \/d 2 - c 2 


t = T - 


Po c 


T = — (t + \/t 2 + c 2 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 


The complex potential function for the horizontal cross flow in the xij-plane is 


W = 3 S \/x 4 2 + £> 2 


where 

P = -^-(d + \jt 2 + d 2 ) 


(38) 


(39) 


Then the conjugate of the complex velocity is obtained by differentiating equa- 
tion (38) by the chain rule and using the transformation equations (29), (31), 
and (32): 
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dW dW dx 4 dx 3 dx 2 

dx-j dx^ dx^ dX 2 dx-j 


3sX4 V X3 Vi 

\|x 4 2 + V*) \Jx 3 2 - d 2 ) \ 


P0 2 \/l 1 X 1 \ 

x 2 2 / \ 2 2 \|x-) 2 - c 2 / 


Consequently, 


v h 


iUh . * v * v, 

2 \£F7vi)\{&T-d*)\ 




(40) 


Resolving the right side of equation (40) analytically into its real and imagi- 
nary parts is a tedious procedure resulting in rather lengthy expressions. 
Calculation of the complex velocity is considerably simplified by retaining 
the complex variables in the programming code for the calculation of dW/dx^ 
and then taking the real and imaginary parts of the calculated value at each 
point . 


CONCLUDING REMARKS 

Some special problems relating to missiles with elliptic-cross-section 
body shapes have been treated by slender body methods. Expressions have been 
derived for the velocity distribution in the nose section when the ellipse 
eccentricity is varying longitudinally on the missile. The cross-flow velocity 
on a triform fin section is also studied. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 29, 1977 
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Figure 1 Nose section geometry with 
longitudinally varying cross-section 
eccentricity. 



Figure 2.- System of axes. 
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Figure 4.- Sequence of contours and parameters 
for mapping missile side flow into flow over 
a flat plate. Numerals denote corresponding 
points . 
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